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Abstract 

Recent experimental measurements of Ag impurity diffusion in the S5 (310) grain 
boundary (GB) in Cu revealed an unusual non-Arrhenius behavior suggestive of a 
possible structural transformation [Divinski et al, Phys. Rev. B 85, 144104 (2012)]. 
On the other hand, atomistic computer simulations have recently discovered phase 

t-h transformations in high-angle GBs in metals [Frolov et al., arXiv:1211.1756v2 (2013)]. 

In this paper we report on atomistic simulations of Ag diffusion and segregation in two 

r — different structural phases of the Cu £5 (310) GB which transform to each other with 

temperature. The obtained excellent agreement with the experimental data validates 
the hypothesis that the unusual diffusion behavior seen in the experiment was caused 
by a phase transformation. The simulations also predict that the low-temperature 
GB phase exhibits a monolayer segregation pattern while the high-temperature phase 
features a bilayer segregation. Together, the simulations and experiment provide the 
first convincing evidence for the existence of structural phase transformations in high- 
angle metallic GBs and demonstrate the possibility of their detection by GB diffusion 

\2 measurements and atomistic simulations. 

Motivation. Structural transformations at grain boundaries (GBs) are of fundamental 
interest as a case of phase transitions in low- dimensional systems. In addition, GB phase 
transformations can have a significant impact on microstructure, mechanical behavior and 
transport properties of polycrystalline materials [H E]. A number of GB phases have been 
found in ceramic materials [3l 0], where they often appear in the form of intergranular 
thin films and are referred to as "complexions" [5]. In metallic alloys, several phases with 
discrete thickness have been observed, such as the segregated bilayer structure believed 



to be responsible for the liquid-layer embrittlement effect [6j. However, despite decades 
of research little is known about GB phase transformations in single-component metals, 
apart from the recently found dislocation pairing transition in low-angle GBs composed of 
discrete dislocations [7] and the prediction of temperature-induced phase transformations 
in high-angle GBs [8]. Presently, there is no direct experimental evidence for structural 
transformations in metallic GBs. The main reason lies in the extreme difficulty of GB 
structure characterization by high-resolution transmission electron microscopy (HRTEM) 
at high temperatures [9j [10] . 

Several indirect methods have been applied to probe the structural state of GBs, includ- 
ing GB diffusion |11J, GB mobility [12] and GB sliding resistance [13]. Recently, impurity 
diffusion of Ag in a Cu bicrystal with a well-characterized £5 (310) GB has been measured 
using a radioactive isotope of Ag (E being the reciprocal density of coincidence sites and 
(310) the GB plane) [Hj. At high temperatures, the GB diffusivity was characterized by 
the diffusion flux P = D g i,sS, where D g b is the GB diffusion coefficient, s is the impu- 
rity segregation factor and 5 the thickness of the GB core considered as a uniform layer. 
The Arrhenius plot shown in Fig. [T] reveals a distinct break in the slope at temperatures 
around 800 to 850 K. This observation correlates with the similar break found previously 
for Au diffusion in a similar GB [15]. Such breaks in the slopes of Arrhenius plots usu- 
ally signify an abrupt change in the diffusion mechanism [11, 16]. On this ground, it was 
suggested [HI [15] that the £5 (310) GB undergoes a structural phase transformation at 
about 800-850 K. While a plausible explanation, its confirmation by direct experimental 
observations is hampered by the mentioned difficulties in HRTEM characterization of GBs 
at high temperatures. 

Much of the current knowledge about GB structures comes from atomistic computer 
simulations [H [2]. Many previous simulations have shown a behavior where GBs exhibit 
structural disorder at high temperatures and eventually melt by either turning into a 
continuously growing liquid film or triggering bulk melting [Tj [2], [HJ [T7J [18]. However, 
until the recent report [8] there was no compelling simulation evidence for transforma- 
tions between different ordered structures in high-angle metallic GBs. Observation of such 
transformations had been precluded by inadequate simulation methodology which prohib- 
ited variations in atomic density in the GB core. When such variations were allowed, 
several alternate structures and first-order phase transitions between them were found in 
the £5 (210) and £5 (310) GBs in FCC metals [8j. For the Cu £5 (310) GB, a transforma- 
tion between two structures was found, with split-kite structural units stable above 800 K 
and kite units stable at lower temperatures (Fig. [2]). The kite structure had been known 
previously, whereas the split-kite structure is new and was first discovered in [8]. 

The finding of this phase transformation [8] and the correlation between its temperature 
and the experimental temperature of the change in the diffusion activation energy [TJ] lend 
more credence to the hypothesis that the observed diffusion behavior was caused by a GB 
phase transformation. However, a convincing proof requires a calculation of the diffusivity 
in impurity-segregated boundaries in the vicinity of the transformation temperature and 
demonstration that it indeed reproduces the experiment. It is the goal of this paper to 
conduct such calculations and validate the proposed hypothesis. We also wish to under- 
stand the effect of the transformation on Ag segregation and the role of this segregation in 



the diffusion behavior. 

Methodology. The molecular dynamics (MD) simulations were performed in the 
canonical (NVT) ensemble and employed the LAMMPS code [19J and the embedded-atom 
potential for copper [20]. The simulation block with a plane £5 (310) GB with the stan- 
dard kite structure was created by the usual geometric construction [21] and contained 
around 47,000 atoms. The block had the approximate dimensions of 6.6x12.5x7.00 nm 3 
with periodic boundary conditions parallel to the GB plane x-z and free surfaces in the 
y direction normal to the GB plane. The split-kite GB structure was first created by an 
isothermal 800 K anneal of a GB terminated at free surfaces. The central part of this block 
was then carved out to create a new simulation block with the same dimensions and the 
same boundary conditions as for the standard kite structure. The obtained GBs (one with 
kites and the other with split kites) were disconnected from sinks and sources of atoms and 
unable to vary their density. Consequently, their structures remained unaltered during the 
subsequent anneals. Using these simulation blocks, diffusivities of the two GB structures 
could be studied separately over the same temperature interval. 

For each GB structure, a series of isothermal anneals was performed at temperatures 
from 750 to 1000 K. The simulation time t increased from 10 ns at 1000 K to 70 ns at 750 
K. Multiple snapshots saved during the simulations were used to analyze the trajectories 
of atoms diffusing within the GB core. The GB self-diffusion coefficient was computed 
from the Einstein relation DL = (x 2 ) jit for the direction parallel to the tilt axis and 
Dg b = (z 2 ) jit for the perpendicular direction. To validate this methodology, we verified 
that the mean-squared displacements indeed increased linearly with time as prescribed by 
the Einstein equation. 

While the obtained self-diffusion coefficients are useful as a guide for future research (see 
discussion below), the experimental measurements [14] were made for Ag impurity diffusion. 
The generalized Fisher model [TT] employed in [14] for the processing and interpretation of 
the experimental concentration curves assumes that the impurity atoms segregate to the GB 
and remain in thermodynamic equilibrium with adjacent lattice regions. It was therefore 
necessary to create equilibrium segregation of Ag in the simulated models. To this end, we 
used the semi-grand canonical Monte Carlo method [22] with Cu-Ag interactions modeled 
with the embedded-atom potential [23J. At each temperature T, the desired equilibrium 
chemical composition c inside the grains was established by adjusting the imposed chemical 
potential difference A/i between Ag and Cu. The lattice constant was chosen to eliminate 
mechanical stresses inside the grains. At each T and A/i, an equilibration anneal was 
performed followed by a production run comprising 4 x 10 5 to 6 x 10 5 Monte Carlo steps 
per atom. Several dilute grain compositions ranging from 0.001 to 0.1 atomic percent of 
Ag were implemented at each temperature. At the temperature of 800 K, a more detailed 
calculation of the segregation isotherm was performed with c varying from 0.001% to 0.6%. 

The blocks with the equilibrium Ag segregation corresponding to the temperatures 
between 750 and 1000 K and the same grain composition c = 0.1% were taken as initial 
configurations for MD simulations of Ag impurity diffusion. As for self-diffusion, D g ^ was 
extracted from mean-squared atomic displacements of Ag atoms in the GB core. Since the 
number of such atoms was significantly less than for self-diffusion, longer MD simulation 
times ranging from 50 to 70 ns were implemented. 



Results. Fig. [3Fa) shows the segregation isotherm at 800 K. The GB concentration c g b 
was obtained by computing the amount of Ag segregation at a constant total number of 
atoms [23] and assuming the GB thickness 5 = 0.5 nm. The segregation factor s = c g b/c 
appearing in the generalized Fisher model [11] must be computed in the dilute regime 
where it is independent of the grain composition. The dilute regime was also implemented 
in the diffusion measurements [14] . where the radio-tracer Ag atoms only served as a probe 
of the GB structure. The dilute segregation factor was obtained by a linear fit through 
the first 3 points near the origin of the plot. Similarly, the segregation factors at all 
other temperatures were calculated using the first 2-3 points where the isotherm was linear 
(Henry-type [DEI]). 

It is interesting to examine the segregation behavior beyond the linear regime. As evi- 
dent from Fig. [3Fa), the segregation becomes non-linear and qualitatively consistent with 
the McLean isotherm [1] between 0.05% and 0.1%. At higher concentrations, however, c g b 
does not saturate as expected from the McLean isotherm but continues to increase. Fur- 
thermore, at about 0.2% the isotherm goes through an inflection point and the increase in 
Cgb accelerates. This behavior is consistent with the experimental isotherms reconstructed 
from GB diffusion data [25] and is usually indicative of multiple types of segregation sites 
and/or segregation to multiple layers. Note that the segregation to the standard kite struc- 
ture is significantly stronger than to the split kites. In addition, the two GB structures 
exhibit different segregation patterns as illustrated in Fig. [3]^b,c). In the case of standard 
kites, the Ag atoms tend to occupy the tip of the kite, whereas in the split kites they segre- 
gate to corners of the kites, creating a bilayer segregation pattern. It should be emphasized 
that neither the monolayer nor the bilayer are completely filled with Ag atoms. In fact, the 
total number of atoms segregated into the bilayer is less than for the monolayer. Based on 
the segregation isotherm one can predict that if the GB structure at a given temperature 
transforms from standard kites to split kites, the equilibrium Ag segregation must abruptly 
drop while the segregation type must change from a monolayer to a bilayer. As already 
mentioned, changes in the segregation type from a monolayer to a bilayer (and possibly 
three or more layers) were observed experimentally in metallic alloys [6]. 

It can also be observed that at c ~ 0.3% the segregation to the split-kite structure 
abruptly increases. The new GB structure which forms above 0.3% is less ordered and 
does not follow the bilayer pattern. This segregation jump is suggestive of a segregation- 
induced structural transformation reminiscent of the one found previously in twist GBs 
[26] . However, this simulation setup does not permit a conclusive identification of the new 
phase because the GB is not free to adjust its total atomic density. A more comprehensive 
study would employ a GB terminated at free surfaces and is beyond the scope of this paper. 

The Arrhenius diagram in Fig. EFa) summarizes the dilute segregation factors for the 
kite and split-kite structures. Both structures accurately follow the Arrhenius law but with 
significantly different slopes. The effective segregation energy extracted from the slopes is 
74 kJ/mole for the kite structure and 40 kJ/mole for the split-kite structure. Included 
are also the experimental segregation factors reported in [13]. It should be noted that 
they were not obtained by direct chemical characterization of the GB but rather back- 
calculated from the diffusivities measured in different kinetic regimes. Nevertheless, given 
that the segregation factors span nearly two orders of magnitude, the agreement between 



the simulations and experiment can be considered favorable. 

Similarly, the diffusion coefficients of Ag in the two GB structures accurately follow the 
Arrhenius law [Fig. Utb)] with the activation energies 149 kJ/mole for the kite structure 
and 93 kJ/mole for the split-kite structure (averaged over the two diffusion directions). The 
diffusion coefficient is lower for the split-kite structure at high temperatures and for the 
standard kite structure at low temperatures. The temperature of the cross-over between the 
two diffusivities correlates well with the experiment (cf. Fig. [I]). Furthermore, extrapolation 
of the kite diffusivity to lower temperatures gives good agreement with the experimental 
D g b values measured in the type-C kinetic regime [14] . 

Putting together the diffusion and segregation data, Fig. [5]shows the Arrhenius diagram 
of the GB flux P = D g i,s5. The first notable feature of this plot is the overall close 
agreement between the experimental measurements [H] and independent simulation data. 
For both GB phases, the computed fluxes P accurately follow the Arrhenius law with 
the activation energies 74 kJ/mole for the kite phase and 52 kJ/mole for the split-kite 
phase (averaged over both directions). For the kite phase stable at low temperatures, this 
activation energy compares reasonably well with the experimental values (67 to 75 kJ/mole 
[14] ) extracted from the low-temperature portion of the experimental Arrhenius diagram. 

The cross-over between the diffusivities P of the two GB phases occurs at the tempera- 
ture of 800 K, which is in excellent agreement with the experiment. This cross-over creates 
the characteristic break in the slope of the Arrhenius plot which is strikingly similar to 
the one reported in [14J. This close agreement strongly suggests that the non- Arrhenius 
diffusion behavior found in the experiment [Hj was caused by the structural phase trans- 
formation in this boundary. 

Conclusions. One general conclusion from this work is that atomistic simulations 
are presently capable of predicting GB diffusion and segregation in good agreement with 
experiment. It has also been demonstrated that our simulation methodology permits dif- 
fusion and segregation calculations for individual GB phases. This opens the door for 
future investigations of the effect of GB phase transformations on thermodynamic (segre- 
gation) and kinetic (diffusion) properties of high-angle GBs in metals by means of atomistic 
modeling. For the particular Cu £5 (310) GB studied here, the simulations predict that 
the structural transformation must cause a break in the activation energy of Ag impurity 
diffusion at about 800 K, with a lower activation energy above this temperature. This pre- 
diction is in excellent agreement with the recent experimental study [14]. This agreement 
confirms that the non-Arrhenius diffusion behavior reported in [T4] indeed constitutes an 
experimental observation of a structural phase transformation in a high-angle GB. 

Our simulations predict that the phase transformation must cause a similar break in 
the activation energy for Cu self-diffusion (Fig. [5]) . This suggests that the non-Arrhenius 
behavior found for Ag diffusion in Cu [H] was primarily caused by the effect of the phase 
transformation on the diffusion coefficient D g b rather than on segregation. Nevertheless, 
segregation also plays an important role as it affects the activation energies of the mea- 
sured diffusion fluxes P in the two phases. Direct segregation measurements may provide 
additional information on the GB transformation. Our study shows that atomistic simula- 
tions are capable of reproducing complex segregation patterns, such as the bilayer, found 
experimentally in metallic alloys [6 J and ceramics [3j H] . This work motivates future mea- 



surements of Cu self-diffusion in the S5 (310) GB to test the predictions shown in Fig. pi 
Cu self-diffusion measurements can employ the radioactive isotopes 64 Cu or 67 Cu and are 
extremely challenging [271 EH] ■ Nevertheless, such measurements could be conducted in the 
near future using bicrystalline samples as it was done in |14j . 
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Figure 1: Arrhenius diagram of impurity diffusion of Ag [14] and Au [15] the Cu E5 (310) 
GB obtained by radio-tracer measurements on bicrystals. For Ag diffusion, the circles and 
squares are for diffusion parallel and normal to the tilt axis, respectively. Note the break 
in the slope of both plots at high temperatures. 
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Figure 2: Phase transformation in the Cu S5 (310) GB at the temperature of 1000 K [8]. 
The initial structure of the boundary was composed of kite-shape structural units (light 
blue). A new phase composed of split-kite structural units (red) grows into the GB and 
eventually penetrates all through the sample. This transformation is enabled by supply of 
atoms from the surface. 
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Figure 3: (a) Isotherm of Ag segregation at the Cu £5 (310) GB at the temperature of 
800 K. The inset shows the low-concentration part used for the calculation of the dilute 
segregation factor. The images show the segregation patterns for the kite (b) and split-kite 
(c) GB structures. For clarity the atomic structures were quenched to K to eliminate the 
thermal noise due to atomic vibrations. 
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Figure 4: Arrhenius diagrams for (a) dilute segregation factor s and (b) diffusion coeffi- 
cient of Ag in the Cu £5 (310) GB. The experimental data [H] is shown for comparison. 
The experimental segregation factors were back-calculated from diffusion data. Diffusion 
parallel and normal to the tilt axis is represented by filled and open symbols, respectively. 
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Figure 5: Arrhenius diagram of the diffusion flux for Ag impurity diffusion (upper curves) 
and Cu self-diffusion (lower curves) in the Cu £5 (310) GB. For Ag diffusion, the experi- 
mental data [14J is shown for comparison. Diffusion parallel and normal to the tilt axis is 
represented by filled and open symbols, respectively 
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